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Condensation Bottleneck Driven by a Hidden Van Hove Singularity as the Origin of 
Pseudogap Physics in Cuprate High-Tc Superconductors 
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We propose a new approach to understand the origin of the pseudogap in the cuprates. The near- 
simultaneous softening of a large number of different ^-bosons yields an ordering bottleneck, wherein 
the growth of magnetic correlations with decreasing temperature is anomalously slow, leading to 
extended ranges of short-range order. This effect is not tied to the Fermi level, but driven by a 
Van Hove singularity (VHS) nesting that is strongest near a temperature Tvhs that scales as the 
energy separation between the Fermi energy and the energy of the VHS peak. By identifying Tvhs 
as the pseudogap onset temperature T*, we explain many characteristic features of the pseudogap, 
including the observed transport anomalies and the termination of the pseudogap phase when Evhs 
crosses the Fermi level. The condensation bottleneck (CB) provides a new pathway for understanding 
strong correlation effects in the cuprates. We find that LSCO lies close to an anomalous disorder-free 
spin glass quantum critical point, where the frustration is due to the CB. To study this point we 
develop an approach to interpolate between different cuprates. 

PACS numbers: 


I. INTRODUCTION 

Evidence is growing that the ‘pseudogap phase’ found in cuprates is actually home to one or more competing phases, 
including a variety of stripe, spin-, or charge-density wave (S/CDW) phases.The CDW phase, in particular, 
has stimulated considerable interestMany of these phases, including superconductivity, seem to appear at 
temperatures well below the pseudogap temperature T*, so the exact relation between the pseudogap and these other 
phases remains elusive. Indeed, the real puzzle is understanding why the pseudogap phase bears so little resemblance 
to a conventional phase transition. Here we demonstrate that the pseudogap phenomenon arises from a condensation 
bottleneck (CB), where a large number of density waves (DWs) with similar g-vectors attempt to soften and condense 
at the same time, leading to anomalously low transition temperatures and extended ranges of short-range order, 
characteristic of the pseudogap phase. We note that many treatments of correlation effects - e.g., path-integral 
approaches - begin by finding a mean-field order parameter at some particular g-vector - a form of random-phase 
approximation (RPA) - and then include the effect of fluctuations. Then if several instabilities compete, the one 
with lowest free energy wins out. In this case, the CB represents a particularly strong correlation effect, insofar as 
one cannot solve the problem one-g at a time, but must account for the g-competition from the start, to account 
for the entropy associated with the g-manifold. We note that RPA-type approximations extend well beyond simple 
mean-field theories, arising in spin-fermion models and in Hertz-Millis^^’^^ type of quantum critical theories when the 
Stoner denominator is approximated by a T-independent Ornstein-Zernicke^^ form, and in all cases they miss the 
essential and strongly T-dependent mode-coupling physics. 

Similar effects have arisen in the past. In Overhauser’s theory of CDWs in alkali metals, the g-vectors for 
all points on a spherical Fermi surface (FS) become unstable simultaneously and cannot be handled one at a time. 
These effects are often referred to as phonon entropy or more generally as boson entropy, and can be analyzed via 
vertex corrections which take proper account of mode coupling^^’^^. In excitonic theories, the various DW instabilities 
represent condensation of some electronic boson at a single g-vector, with triplet excitons corresponding to SDWs and 
singlet excitons to CDWs^^. Summing only bubble and ladder diagrams reproduces the RPA results with BCS-like 
ratios of gap to critical temperature, 2 A//cbTc- Going beyond ladders (e.g., Bethe-Salpeter equation) incorporates 
mode coupling with enhanced A/Tc ratios (compare Refs. [31] and [32]; for a review of excitonic insulators, see 
Ref. [30]). A particularly interesting analogy is provided by the Bose condensation of excitons. The mean-field 
transition is found to correspond to the temperature at which excitons are formed. However, since the excitons are 
localized in real space, they are greatly spread in g. When fluctuations are included, the real transition lies at a much 
lower T, when all excitons condense into the lowest g state. Here we develop a similar theory for the cuprates via a 
self-consistent renormalization calculation of the vertex corrections. We find that the CB not only drives pseudogap 
physics, but can potentially lead to a novel disorder-free spin-glass phase. Similar effects are likely to be present in 
many other families of correlated materials. 

We find that the phase diagram of La 2 -ccSra,Cu 04 (LSCO) is different from that of other cuprates. To explore 
the transition between them, we simplify the first principles dispersions to equivalent ‘reference family’ dispersions, 
depending only on the three nearest neighbor hopping parameters, t, t', and t" (Supplementary Material Section I). 
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We find that all cuprates fall along a single cut in t' jt — t'Yt-space, namely t" — —t'{2, with LSCO characterized by 
a smaller value of t'. 


II. RESULTS 

A. Bosonic VHS and Origin of the Susceptibility Plateau 

A proper understanding of the CB involves two different issues: firstly, since Fermi surface (FS) nesting typically 
singles out only a few discrete g-vectors, what causes many g-vectors to soften together? Secondly, a formalism needs 
to be developed to handle the effects of bottleneck. Both of these issues are addressed here. Remarkably, in cuprates 
we find that the underlying source of strong mode coupling can be traced back to the Lindhard susceptibility xo- 
In a typical calculation of classical or quantum phase transitions^^’^^, material parameters are introduced via xo, 
which is used to define an interacting susceptibility x- In onr calculation, xo is calculated from density functional 
theory (DFT) bands corrected by a GW self-energy [Methods Section], and x is the resulting RPA susceptibility. The 
phase transition then corresponds to the vanishing of the denominator of x at frequency cj = 0 for some momentum q. 
To properly incorporate fluctuation effects, this ‘Stoner denominator’ is then typically reduced to Ornstein-Zernicke 
(OZ)^^ form. 
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in terms of various deviations from the critical point (in g, w, and a ‘tuning parameter’ which is proportional to 
where ^ is the correlation length). Here z is a. dynamic exponent and q is measured from the Q where the susceptibility 
has a peak. An important result of our calculation is that in cuprates, this OZ form must be treated carefully, since 
the coefficients of the deviation parameters can be strong functions of temperature and doping. In particular, we 
find that the susceptibility inverse curvature (coefficient of g^) can diverge due to a competition between conventional 
Fermi-surface nesting and Van Hove singularity (VHS) nesting. 

Whereas many properties of a Landau Fermi liquid are determined by their values at the Fermi level, the suscep¬ 
tibility is an exception, having both bulk and Fermi surface contributions. While the FS part contributes a ridge 
to Xo that is a map of the FS at g = 2kF^'^ (where /ci? is the Fermi wave vector), in the cuprates there is also an 
important bulk contribution, which provides a smoothly varying background, peaking at (tt, tt) and giving rise to the 
near-(7r, 7r)-plateau in the susceptibility. This peak can shift the balance of the FS nesting to g-vectors closer to the 
peaks, and in special cases can lead to commensurate nesting away from the FS nesting vector, at exactly (7r,7r). 
Moreover, as T increases, coherent FS features are washed out, leaving behind only the commensurate bulk contribu¬ 
tion. This peak is a bosonic VHS (b-VHS), the finite-g analog of the Van Hove excitons found in optical spectra^^, 
but present in the intraband susceptibility. However, it is a “hidden” b-VHS. Despite the fact that it is pinned to 
zero energy independent of doping or hopping parameters^ it is hidden in the sense that the effective density-of-states 
(DOS) exactly at the b-VHS peak almost always vanishes. Some consequences of this are discussed in Supplementary 
Materials Section H. 

The imaginary part of the susceptibility can be thought of as the DOS of electronic bosons, electron-hole (e-h) 
pairs, which may become excitons when a Coulomb interaction is turned on. If the renormalized dispersion of a single 
electron is Ck with wave vector k, then an e-h pair at wave vector g has a dispersion ujq{k) = ek-\-q — = ~2eg_(/c), 

where eq±{k) = (e/^ ±e/c+g)/2, and a Pauli blocking factor Afj^^q = /(e/^+g) — /(e/^). Then the corresponding pair DOS 
is Dq{uJ) = Afk^qS{uj — O 0 q{k)) = Xo{Qi^)' Foi* LSCO, the dominant pairs are those at g = Q = (7r,7r), the AF 
nesting vector. The associated dispersion co’g(7r,7r), plotted in Fig. 1(a), resembles the electronic dispersion e/^, but 
with an important distinction: it depends only on eg_(/c), whereas all of the hopping terms that shift the electronic 
VHS away from half filling (P, t") are contained in eg+(/c), i.e. the b-VHS is pinned at cj = 0. Since Xo is an odd 
function of cj, XoiQi^ = 0) = 0. However, excitonic instabilities depend on a Stoner criterion, and hence on 


Xo(9,w 
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The FS contribution thus arises from uo' ^ 0, while bulk contributions arise from peaks in the bosonic DOS, such 
as the b-VHS peak in Xo near (7r,7r). However, while the b-VHS is pinned at cj = 0, Fig. 1(a), its weight vanishes 
at T = 0, due to the Pauli blocking factor, Af^^g = 0 near k = (7r,0) at T = 0, Fig. 1(b). Finite T restores weight, 
optimally near lOOOK, although A/ always vanishes exactly at (7r,0). In the Supplementary Materials, Section H, we 
deconvolve the near-(7r,7r) susceptibility to show that it is a superposition of bulk and Fermi surface features. The 
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FIG. 1: Coherent-incoherent crossover at the VHS. (a) Pair dispersion, uq, as a function of k, and (b) the corresponding 
weight A/, at T = OK (red), 500K (blue), and 2000K (violet) for LSCO ( x = 0). (c) Phase diagram for LSCO, showing 
VHS-related crossovers and regions where the susceptibility peaks at different ^-vectors. Here xq is the doping at T = 0 
(plots are at constant Ep), but for qualitative purposes we can assume x xq. Dominant ffuctuations are at commensurate 
(7r,7r) (white shaded region) or incommensurate (7r,7r — (5) (green shaded region); crossovers are Tvhs (red short-dashed line), 
DOS peak (green long-dashed line), coherent-incoherent crossover (pink shaded region), (violet dot-dashed curve), and the 
position of the ( 7 r, 7 r) peak vs doping (light blue dot-dashed line) - the approximate electron-hole symmetry point. 


bulk feature is dominated by the b-VHS, while the FS feature smears out with increasing T in a coherent-incoherent 
crossover, pink shaded region in Fig. 1 (c). 

To demonstrate the close similarity to VHS effects, we also plot two characteristic features of the VHS. The crossover 
scales with both ksTyns = Ep — Eynsi dark red dashed curve in Fig. 1 (c), and with Ty (violet dot-dashed line), the 
temperature at which the Sommerfeld constant 7 = dSjdT^ has a peak. Since the Fermi function evolves smoothly 
with T, we wxpect sharp Fermi surface features to wash out as T is increased, but why should this coherent-incoherent 
crossover move to lower T as doping is increased? Because a source of entropy, the b-VHS, is moving closer to Ep- To 
demonstrate the connection with entropy, we look at the Sommerfeld constant 7. We calculate 7 from the electronic 
dispersion, assumng a paramagnetic phase, to avoid compications arising from phase transitions. At T = 0 , 7 is 
proportional to the DOS, and hence diverges when the VHS crosses Ep^ At finite T the peak in 7 represents the 
excess entropy associated with mode coupling. The other curves in Fig. 1 (c) will be considered in the Discussion 
Section below. 


B. Strong Mode Coupling Leads to an Extended Range of Short-range Order 

In conventional mode-coupling calculations^^’^^'^^, the OZ parameters are assumed to be T- and doping- 
independent, and the resulting physics becomes quite simple. For 2D materials, the Mermin-Wagner (MW) theorem^^ 
is satisfied, and the mean-field transition at turns into a pseudogap onset at T* ^ with a crossover to long- 
range order when interlayer coupling is strong enough - in short, not much changes from the mean-field results. 
However, when we incorporate realistic susceptibilities we find that this OZ form fails to properly account for the 
strong mode-coupling effects, leading to dramatically different results. In particular, the entropic effects are encoded 
in a strong T-dependence of the model parameters, which greatly slows down the growth of correlations. This in 
turn leads to an extended T-range of short-range fluctuations, which is typical of pseudogap physics. We find that, 
particularly at high T, the doping dependence reflects the evolution of the anomalous VH scattering. We note that 
while the OZ form is used in quantum critical theory^^’^"^, it is explicitly stated that it is to be used only in a limited 
T-range, and only in the absence of FS nesting^^, both of which are violated here. 

Mode coupling modifies the bare Coulomb interaction V, producing a vertex-renormalized effective V, Ugp = TU 
with F = 1/(1 + A), while A is found self-consistently from [see Methods section below] 

A = (3) 

JUc A_ — Usp 

where Aq = 12i4/xo(Qo5 0): Xo(Qo5 0) = 0)]^ is the mode-coupling parameter, and Uc = l/xo(Qo5 0)- 

Here, Xo^ ^ denoted by V_, is the variable of integration and we introduce a corresponding susceptibility density of 
states (SDOS) V_. In Section III of Supplementary Materials, we show that Eq. 3 is closely related to excitonic Bose 
condensation, with A proportional to the effective number of bosons. 
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FIG. 2: Structures in LSCO correlation length, (a) SDOS for undoped LSCO at several temperatures; light-blue dotted 
line gives OZ form of SDOS. (b) Corresponding temperature dependence of correlation length ^th- For x = 0, the red solid line 
with hlled dots is ^th calculated from the solid lines in (a), while the thin green dot-dashed line is based on the dashed lines in 
(a), showing that ^th is relatively insensitive to the structures away from threshold. In contrast, the blue dashed line is based 
on the shape of the SDOS at T = 2000i^", but shifted and renormalized to match the SDOS at lower T, illustrating sensitivity 
to the leading edge structure. The black dotted line illustrates the scaling ^th oc (c) Calculated ^th replotted for x = 0 

(red solid line with hlled circles) compared with ^ (blue hlled circles) and with experiment (green dot-dot-dashed hne)^°; the 
blue dotted line is twice ^th- 


Figure 2 illustrates the profound effects that strong mode coupling has in LSCO, as well as the complete inability 
of the OZ approximation to capture this physics. The SDOS, Fig. 2(a), contains VHS-like features characteristic of 
conventional DOSs. However, the singular behavior of Eq. 3 involves only features near threshold, Uc = min{X_), 
which evolve strongly with T, see inset to Fig. 2(a). For T > 0, the threshold behavior is always a step at Uc, 
indicative of a parabolic peak in xo with curvature inversely proportional to the step height. Thus for a qualitative 
understanding of the evolution of ^ (Eq. 1) with T we can assume an OZ form of x, but with a strongly T-dependent 
step height A 2 {T). This allows a threshold correlation length (red solid lines in Eigs. 2(b) and (c)) to be defined 
from Eq. 1 as the inverse half-width in q, assuming that the curvature (oc is independent. In reality, the 

curvature increases with g, and the correlation length ^ obtained from the renormalized susceptibility half-width is 
typically a factor of 2 larger. Eig. 2(c) compares the x = 0 values of ^th (solid red line) and ^ (filled blue circles) 
vs 1/T; the values of ^ lie on the blue dotted curve, which represents twice and are in good agreement with 
experiment (green dot-dot-dashed line)^^. 

Eigure 2(b) illustrates how strong mode coupling slows down the correlation length divergence. Undoped LSCO 
(red solid line) shows two regions I and III of exponential growth of ^th with decreasing T, separated by an anomalous 
region II where ^th actually decreases with decreasing T. We will not discuss the MW-like divergence at low T (region 
III). In the high-T limit (region I), the leading-edge parabolic curvature is quite small, and if it were T-independent, as 
in the OZ approximation, the growth in would follow the blue dashed line (Eq. 11 below), but thermal broadening 
causes the curvature to decrease with increasing T, leading to the faster growth of the red solid line. The anomalous 
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FIG. 3: t’ dependence of correlation bottleneck, (a) ^th for the minimal reference family at several values of t'/t [see 
legend] and t" = 0. (b-d): SDOS at three values of t'/t = -0.21 (b) [as in Fig. 2(b)], -0.27 (c), and -0.35 (d). All curves are 
shifted to line up the (7r,7r) data at AX- = X- — A_(7r, tt) = 0. (e) ^th for the PA reference family at several values of t'/t 
[see legend] and t" — —t'l2. 


behavior in region II will be discussed next. 


C. Exploring parameter space 

To understand the origin of the CB in region II, it is necessary to explore hopping parameter space away from the 
physical cuprates. To do this, we introduce the notion of reference states: states with simplified hopping parameters 
(only t, t', and t" nonzero) but which match the phase diagram of the real cuprates in a well-defined way [Supplemen¬ 
tary Material Section I]. This allows us to tune the system between LSCO and Bi2212, and explore phase space beyond 
these limits. We study two important cuts in f /t — t"/t space, a minimal cut it" = 0) and the Pavarini-Andersen 
[PA] cut it" = —t'12) - the latter seems to best capture the physics of the cuprates. By tuning t' we unveil the origin 
of the CB as a localization-delocalization crossover tied to the crossover from (7r,7r)- to FS-nesting. As a byproduct, 
we gain insight into why LSCO is so different from other cuprates, and how cuprates evolve from the pure Hubbard 
limit (t' = 0). 

Figure 3(a) shows the T-evolution of ^th for several values of t' along the minimal cut, including the data of Fig. 2(b). 
For t'/t > —0.17, the system is characterized by commensurate (C) (7r,7r) order with a finite Neel temperature 
Tat ^lOOOK [the correlation length grows so rapidly that interlayer correlations will drive a transition to full 3D 
order]. Similarly, for t'/t < —0.345 there is incommensurate (I) (7r,7r — (i) order with Tyr about a factor of 10 smaller. 
But for intermediate t'/t the C-I transition is highly anomalous, with correlation length orders of magnitude smaller 
than expected. Figure 3(b) shows that a similar evolution follows along the PA cut in parameter space. The reason 
for this anomalous behavior can be seen by looking at the leading edge SDOS in the crossover regime, shown for three 
values of t' along the minimal cut in Figs. 3(c)-(e). For ease in viewing, these curves have been shifted to line up 
the SDOS at (tt, tt) at all T. It is seen that the anomalous collapse of ^th is associated with a rapid growth of the 
step height A 2 , culminating in a near-divergence at = —0.27t, where the leading edge curvature goes to zero. This 
divergence coincides with the C-I crossover of the leading edge SDOS [Supplementary Materials Section I]. Here many 
different g-vectors compete simultaneously, frustrating the divergence of any particular mode. This is the electronic 
analog of McMillan’s phonon entropy: if many phonons are simultaneously excited, the transition is suppressed to 
anomalously low temperatures. Note that ^th drops by 9 orders of magnitude at T = 200iF when t'/t changes from 
-0.17 to -0.27, then grows by a similar amount at lOOK when t'/t changes from -0.27 to -0.345. [Over this same 
range, the true ^ will be frozen at the value corresponding to the half-width of the (tt, tt) plateau.] The green shaded 
region in Fig. 3(c) shows that the range of the anomalous growth H of ^ in Fig. 3(a) coincides with the range of rapid 
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growth of the SDOS leading edge; such behavior is absent if a T-independent OZ form (black dashed line) is assumed. 
The divergence of the SDOS at T = 0 signals that the system at t' = —0.27t is highly anomalous, with an infinitely 
degenerate ground state. This represents an anomalous form of spin glass arising in the absence of disorder, with the 
frustration arising from strong mode coupling. Note also that when the C-I transition is at T > 300K, ^(T) has a 
sharp downward cusp at the transition, while above the transition ^ < a is strongly suppressed. 

From the relationship between ^th and the step height A 2 , Eq. 11 below, we see that a diverging A 2 will cause 
^th 0- The divergence arises when the susceptibility at (tt, tt) crosses over from a maximum to a local minimum. 
In the latter case, the maximum intensity of x(g, 0) is spread along a ‘ring’ in g-space surrounding (7r,7r), so that 
the inverse susceptibility resembles a ‘Mexican hat’. Approximating the x“^-dispersion by a Mexican hat form 
— A 2 q^ + A/^q^ leads to a threshold divergence N- (x ^ — Uc) Fig. 2(f). Note that the SDOS-divergence 
resembles a ID VHS in the conventional DOS, even though here the 1-D direction is the radial direction away from 
(7r,7r). The divergence of x along a ring can be thought of as a 2D analog of Overhauser’s effect^^. Similar effects 
are found in Bi2201 [Supplementary Materials Section IV]. Finally, for t' near but greater than the divergence is 
avoided, but proximity to the C-I transition leads to the rapid growth of a narrow peak in the leading-edge SDOS, 
Fig. 2(e). This in turn can cause an approximate power-law growth of ^th ^ (black dotted line in Fig. 2(b)), 

as calculated in the Methods section below, although the sharp turn-on of structure in the SDOS actually leads to an 
anomalous decrease in with decreasing T. 


III. DISCUSSION 
A. Cuprate Pseudogap 

Figure 2(g) captures an essential aspect of pseudogap physics: an extended regime of phase space where correlations 
remain only short ranged. The reason that the C-I transition is so anomalous is that it is also an ineoherent-eoherent 
transition, with the bare susceptibility for small \t'\ dominated, not by FS nesting, but by a broad peak at (7r,7r) 
associated with a b-VHS, crossing over to conventional (coherent) FS nesting only for t' < t'^ or for larger x. Mode 
coupling associated with self-consistent vertex corrections leads to a regime where the local and extended behaviors 
become strongly entangled, leading to a collapse of the magnetic correlation length and very slow growth of magnetic 
fluctuations. 

While the CB was explained by tuning parameter space, the results bear a striking resemblance to the temperature 
and doping evolution of the pseudogap phase. Experimentally, both LSCO and YBCO display C (tt, tt) long-range 
AFM order at x = 0, with dropping rapidly with doping below x = 0.03, crossing over into a regime of I (tt, tt — ^) 
fluctuations and low-T spin glass behavior. Given the proximity of LSCO to the C-I transition, finding a quantitative 
model of its doping phase diagram may prove difficult. We can however note one plausible scenario. It has been 
predicted^^ that a large Hubbard U will renormalize t' to smaller values at half filling, although the magnitude of the 
effect is debated. If doping screens U, and restores a larger |C|, it could drive LSCO across the anomalous regime II, 
causing T/v to drop by an order of magnitude. Moreover, Fig. 2(g) shows that the C-I transition involves a highly 
disordered regime separating two well-ordered phases. Since ordering tends to lower the free energy of the electronic 
system, the disordered regime may represent a state of high free energy, and doping across this region can lead to a 
regime of [nanoscale] phase separation (NFS), which in LSCO is manifest as the stripe phase. Hence, a model of the 
C-I transition, incorporating effects of NFS and disorder, could provide a good description of LSCO. 

This broad T- and doping-regime of only short-range SDW order is a key characteristic of pseudogap physics. 
Indeed, Tyns is typically close to the measured pseudogap temperature T* in most cuprates. Identifying T* with 
Tvhs can explain a number of puzzling features. In particular, it has been found that the pseudogaps in Bi220I and 
Bi22I2 terminate when Tyns 0? i-c., at the conventional VHS,^^'^^ similar to the crossovers seen in Fig. 1(c), as 
well as in Fig. I of the Supplementary Materials. 

The presence of an excitonic VHS contribution to the incoherent susceptibility has a number of further consequences 
for cuprate physics. First, the VHS onset near Tyns combined with the C-I crossover at Tcoh ^ TyHs/^ can explain 
the anomalous transport properties found near the pseudogap. Thus, the pink shaded region in Fig. 1(c) shows the 
crossover from incoherent susceptibility dominated by the VHS at high-T to coherent, FS-dominated susceptibility 
at low-T. This parallels transport, where for T > T*, the resistivity p varies linearly with Tfi behavior expected 
near a VHS^^. For lower T, the resistivity is mixed, but p ^ T‘^ ^ as expected for a coherent Fermi liquid, is found 
below a Tcoh < Tyns- This picture bears a resemblance to the Barzykin-Fines model of the cuprate pseudogap, 
identifying Tyns and Tyns/^ with T* and T*/3 in their model. The underlying physics of the incoherent-to-coherent 
crossover in their model is related to Kondo lattice physics^^, with the VHS peak standing in for the Kondo resonance 
(see Supplementary Materials Section V). This raises the question of whether a similar mode-coupling calculation in 
heavy-fermion compounds could lead to a similar anomalous entanglement at the f-electron C-I transition. 
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The excitonic instability should be maximal when the susceptibility is approximately electron-hole (e-h) symmetrical 
in doping, and at T = 0 this happens when the VHS is at the Fermi level.At finite T, the e-h symmetrical point 
continues to coincide with a peak in the (tt, tt) susceptibility at doping Xp, but Xp shows a remarkably rapid evolution 
with T towards half-filling (light-blue dot-dashed line in Fig. 1(c)), nearly coinciding with the coherent-incoherent 
crossover. We interpret this temperature dependence as follows: The parameters t' and t" are relevant perturbations 
shifting the VHS from the pure Hubbard model value at x = 0 where t' = t" = 0. When ksT > |t'|, these 
perturbations become irrelevant, so that the effective VHS becomes electron-hole symmetric at x = 0. Lastly, the 
lower branch of the (tt, tt) — (tt, tt —^) commensurate-incommensurate transition (green shaded region in Fig. 1(c)) also 
scales, following Tyns very closely in this simple (t — t' only) model. Indeed, at T = 0, this transition corresponds 
to the Fermi energy falling off of the (tt, tt) plateau, providing another indication that the effective VHS is shifting 
towards half-filling as T increases. In Supplemental Material Section H we show that this is another consequence of 
Pauli blocking. 

The rapid thermal evolution of the (tt, 7r)-VHS peak should be contrasted with the much smaller change in conven¬ 
tional VHS effects found near F. In particular, the DOS peak is dominated by near-FS physics, and hence displays a 
much weaker doping dependence, green long-dashed line in Fig. 1(c). 

An additional consequence is that in cuprates with small t' hopping, such as LSCO, the anomalous VHS suscep¬ 
tibility can dominate even at T = 0, leading to strong deviations from FS-nesting. Thus, in most of Fig. 1(c), the 
peak susceptibility is associated with VHS nesting at q = (7r,7r), and only in an intermediate T-regime is FS nest¬ 
ing is found at — S) (green shaded region). The resulting VHS-FS nesting competition plays a strong role in 
underdoped LSCO, associated with retrograde correlation length change with T, region H of Fig. 2(b), which is a 
signal of proximity to a novel disorder-free spin-glass QCP. Notably, in LSCO commensurate (tt, tt) order disappears 
rapidly by ^ 2% doping, being replaced by incommensurate magnetic fluctuations and low-T spin-glass effects, while 
neutron scattering has found that doped LSCO is close to a magnetic QCP.^^ The very different situation in most 
other cuprates is discussed in Sections I, IV of Supplementary Materials. 


B. Strong coupling physics 

We recall that our self-energy formalism^^ is able to reproduce most spectral features of the insulating cuprates 
in terms of a (7r,7r) ordered phase. We reproduce not only the photoemission dispersions, limited to the lower 
Hubbard bands, but optical and x-ray spectra which depend sensitively on the Mott gap. In a related 3-band model, 
we reproduced the Zhang-Rice result that the first doped holes are predominantly oxygen character, and our overall 
dispersions at half filling agree with [subsequent] DMFT results at least as well as DMFT results from different groups 
agree. [See further Supplementary Materials VL] The problem with our earlier calculations is that they predict long- 
range (tt, tt) AFM order at too high T. The new self-consistent renormalization calculations have only short range 
order, in which case the upper and lower Hubbard band dispersion is reproduced, with a broadening r\j 

An antiferromagnet is considered to be weakly coupled if the Neel temperature increases with V, and strongly 
coupled if T/v decreases with increasing U [T/v ^ J = 4t^/I/]. Since this is a finite T criterion, it is sensitive to the 
boson entropy effects we have been discussing. Indeed, for large U the system is nearly localized and the physics 
is again reminiscent of Bose condensation. The simultaneous softening of many magnetic modes can be readily 
demonstrated from a simple Hartree-Fock model of an antiferromagnet. For arbitrary g, the resulting gap for large U 
is 


/^E= yfm + el - U + Je'i, (4) 

where e_ = (e/^ — ek-\-q)l2 = —2te_. Thus, the difference in energy between any two spin configurations is a quantity 
of order J, so that when T ^ J, the system gains entropy by mixing different g-states, and long-range AF order is 
destroyed.While this effect can be recognized in HF, only a theory that properly accounts for the mode competition 
can resolve it. The present mode coupling model properly accounts for this effect via the vertex renormalization 
factor A, which is an integral over the RPA y. This can be seen from Fig. 2(c), where the measured ^ in LSCO 
(green dot-dot-dashed line)^^ is compared to our calculation (filled blue dots). The experimental data are shown only 
above 300K, since at lower T interlayer coupling drives a transition to long-range order. In the Heisenberg model, 
^ is a function only of T/ J, and these data have been used to measure the exchange J. Hence, our mode coupling 
calculation successfully reproduces this strong coupling result, as was found earlier for electron-doped cuprates^^. The 
above results suggest that the C-I, incoherent-coherent transition is also a localization-delocalization transition, with 
the mode coupling associated with small effective hopping. 

For more insight into the strong coupling limit, we note that in the two-particle self-consistent approach^^, A is 
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determined by a sum rule involving double occupancy, which is fixed by assuming 

Uap/U =< > /(n/2)2, 

which leads to a saturation of Ugp as ^ oo (or < 1/^)- calculation, this saturation arises naturally, 

since Usp can never exceed Uc- 

Quantum critical points (QCPs) in strongly correlated materials are often discussed in terms of deconfined QCPs 
(DQCPs), involving a non-Landau transition between two types of competing order, where a new form of excitation 
emerges exactly at the DQCP. This has been refined for cuprates into an underlying competition between Mott physics 
and Fermi liquid physics, masked by a low-energy order parameter of the FL^^. This is an apt description of the 
current results, with the Mott physics evolving into a low-T AFM and the FL to a spin-density wave (SDW), and a 
spin glass at the DQPT. This confirms the finding from dynamical cluster approximation calculations of the strong 
role of the VHS in Mott physics. 


C. Conclusions 

In conclusion, we find that the pseudogap is driven by tendencies toward magnetic order at or near q = (7r,7r). 
Strong mode coupling effects drive the characteristic anomalies of the pseudogap, in particular the existence of a broad 
doping- and temperature-range of only short- to intermediate-range order. Finally, the underlying cause of the strong 
mode coupling is localization in particular the fact that the extended-to-localized transition is not a simple crossover 
but a competition between FS-dominated and non-FS dominated (VHS-dominated) physics. This competition across 
a manifold of g-states leads to a condensation bottlneck, which requires a new formalism for dealing with multi-g 
mode softening. The close similarity to heavy Fermion physics should be noted. 

This can be restated slightly differently. The original (t-only) Hubbard model is particularly difficult to solve, 
as three separate instabilities are simultaneously present at half filling: the Mott instability, FS nesting, and VHS 
nesting. A finite t' shifts the latter two to finite doping. For large |t'| the three instabilities become well separated, 
but for small \t'\ the two nesting instabilities overlap and compete, leading to frustration and pseudogap physics. 
Significantly, we approach the problem from the intermediate-coupling side, suggesting that the full crossover from 
half-filling to large doping could be explored. 

In turn, this sheds light on strong correlation effects in the cuprates, in particular Mott vs Slater physics. The 
localization associated with Mott physics has several manifestations. One is that coherent hopping is restricted to 
shorter range. That is consistent with stronger effects found at smaller t'. Particularly for Mott physics, a second 
aspect is the breakdown of hybridization. The Cu-0 hybridization in the cuprates spoils Mott physics [half filling does 
not imply one electron on each Cu]. Zhang-Rice singlet formation can be considered as a breakdown of hybridization 
[pure oxygen states at the top of the lower magnetic band, pure Cu at the bottom of the upper magnetic band], 
and is found more generally for (7r,7r) AFM fluctuations near half filling in a three-band model^^. Hence the strong 
evolution with t' from localized to extended physics, with only LSCO close to Mott physics. 

Rice et recently noted: “The need for theoretical methods to handle ... short range correlations ... is a key 
challenge for the future.” Along this line, the present approach appears to explain a number of anomalous features 
associated with the pseudogap, including a coherent crossover associated with a competition between two DW orders 
(see Supplementary Materials Section IV). In fact, we have incorporated an important ingredient for strong-coupling 
calculations: a model of short-range AF order from which a t — J model can be derived. 

In Anderson’s RVB picture, he envisaged a regime of Mott physics where the FS played a neglible role; instead, 
most experimental studies find clear evidence for a well-defined FS, with competing phases associated with FS nesting, 
and a QCP associated with FS reconstruction. Our results have a strong bearing on these competing scenarios: in 
parameter spaee^ t' is a relevant parameter tuning the system away from quasilocalized physics near t' = 0, where the 
FS has a negligible role, into a delocalized regime dominated by FS physics. The crossover is marked by a regime of 
strong disorder, where the correlation length remains small down to very low T. Most families of cuprates lie on the 
delocalized side of this crossover, and hence are most easily understood from a FL-type picture. LSCO appears to 
be on the localized side of the crossover, but close to it, consistent with experiment.We note that the disorder line 
may be hard to directly access: since ordering transitions tend to lower the free energy, the line will be a high-energy 
state with low energy states to either side of it, suggesting an instability to nanoscale phase separation. 

The VHS has been predicted to play a significant role in many materials, particularly in lower dimensional systems 
where y diverges, but clear evidence for this remains sparse. Thus, VH nesting was introduced as a possible cause 
of CDWs in dichalcogenides^^, but this interpretation remains disputed. Even in VO 2 , the striking metal-insulator 
transition has been found to be driven by large phonon entropy.Clearly, one problem is that the VHS has been 
assumed to play a role only when it is near the Fermi level, whereas Fig. 1(c) demonstrates that its influence extends 
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over a much wider doping range. For those who think condensed matter physics lies on the surface of the Fermi sea, 
the VHS is the iceberg, lurking. 

It will be interesting to examine how the present results are modified by effects of disorder and interlayer coupling. 
An important issue is the extent to which the total number of QPs is conserved. More specifically, what is the relation 
between the number of incoherent electrons (1 — Z)N and the effective number of excitons N^ffl Proximity to a 
spin glass phase could explain the failure of recent attempts to create an ‘artificial cuprate from related nickelate 
compounds, and may suggest a path to improved analog states. 


IV. METHODS 


Our calculation is a form of many-body perturbation theory (MBPT) based on Hedin’s scheme. The scheme 
involves four elements: electrons are described by Green’s functions G with DFT-based dispersions renormalized by a 
self-energy E; electronic bosons [electron-hole pairs] are described by a spectral weight [susceptibility] y renormalized 
by vertex corrections F. Neglecting vertex corrections, the self energy can be calculated as a convolution of G and 
W = the GW approximation. This approach has been used to solve the energy gap problem in semiconductors, 
where the F correction leads to excitons via the solution of a Bethe-Salpeter equation, and in extending DMFT 
calculations to incorporate more correlations [e.g., DMFT+GW, etc.]. Our approach here is to extend our previous 
GW calculations [quasiparticle-GW or QPGW^^] to include vertex corrections. 

In QPGW we introduce an auxiliary function Gz = Zj{uj — where the dressed, or QP dispersion is = 

Z^dft^ and is the bare, or DFT dispersion. Gz behaves like the Green’s function of a Landau-type QP - 

a free electron with renormalized parameters that describes the low-energy dressed electronic excitations. However, 
this is a non-Fermi liquid type QP, since the frequency-integral of Im{Gz) is Z and not unity. That is, the Z-QP 
describes only the coherent part of the electronic dispersion, and is not in a 1:1 correspondence with the original 
electrons. The importance of such a correction can be readily demonstrated. Since a Z-QP has only the weight Z 
of a regular electron, the susceptibility [a convolution of two Gs] is weaker by a factor of Z^ than an ordinary bare 
susceptibility. To match this effect in the Stoner criterion requires introducing an effective Ueff = ZU. In contrast, 
MBPT calculations in semiconductors typically set the GW-corrected Green’s function to G^^ = cj —where 
is the average GW-renormalized dispersion^^, thereby missing the reduced spectral weight of the low-energy, coherent 
part of the band. 

We work in a purely magnetic sector, with a single Hubbard U controlling all fluctuations; there is a competition 
in this model between near-nodal (NNN) and antinodal nesting (ANN) which mimics the SDW-CDW competition in 
cuprates, sharing the same nesting vectors^^’^^. 

The self-consistent parameter A is found from a Matsubara sum of the susceptibility 


where N is the number of g-points. 


A = 


AoT 

N 




x(q, i^n) 


Xo(q,iu)n) 

1 + A - C/xo(q,'jw„)’ 


( 5 ) 


(6) 


and the summation in Eq. 5 can be transformed: 

T . . . f cPqa‘^ duj cj 


N 


/ • A f ^ J.7 / ^ A /// • TA 

2^ X(q,^^n) = J 


with a the in-plane lattice constant. 


— + ^2^, 


^ f (Pqa^ /*“ du „ 


- f cPqa? ,, . 

The term Ai introduces a small, nonsingular correction^^ to Eq. 5, which we neglect. 


( 7 ) 

(8) 

(9) 
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FIG. 4: Origin of SDOS features. ( a) Inverse susceptibility Xo ^(q) Bi2201 at doping x = 0. Note that Uc = mm(xo 

(circled region). Curves are at T = 10 (red) and lOOK (blue), (b) Corresponding susceptibility density of states (SDOS) 
N-{X-)^ plotted horizontally, corresponding to the susceptibility of (a). (c,d) Blowups of circled regions in (a,b). 


In order to explore the role of susceptibility plateaus for realistic band dispersions, the self-consistency equation is 
evaluated numerically. For this purpose, we note that since x^^ has dimensions of energy, a plot of 0) resembles 

a dispersion map. Hence we can define a SDOS: 

j^ = j N-{X-)dX.. (10) 

Figure 4 illustrates how integrating over the inverse susceptibility. Fig. 4(a), leads to the SDOS, Fig. 4(b), at x = 
0 where the NNN plateau is dominant. Here we use hopping parameters appropriate to the DFT dispersion of 
Bi2Sr2Cu06+a: (Bi2201), t = 419.5, t' = —108.2 and t" = 54.1 eV (Supplementary Materials Section IV). This yields 
a doping phase diagram that is qualitatively similar to that of most cuprates, except LSCO. By comparing Figs 4(a) 
and 4(b), one can see how features in x^^ translate into features in V_. Thus, the intense, flat-topped peak in N_ at 
small values of x^^ represents the NNN plateau. Its broad leading edge (smaller x^^) is controlled by anisotropy of 
the plateau edge between (tt, t: — 5) and (tt — (5, tt — (5), while the sharp trailing edge corresponds to the local maximum 
of Xo ^ hor the ANN peak, its leading edge scarcely leaves any feature in V-, but a local maximum translates 

into a large peak in V_. Finite temperature, T = lOOK (blue line) rounds off the cusp in hig. 4(c), leading to a 
step in V_, Fig. 4(d), but otherwise has little effect. 

Once N_ has been calculated, Eq. 3 can be evaluated numerically. The singular part of the integral is treated 
analytically, and the remainder numerically, with u and U approximated as constants, u = 0.8eV“^, while U = 2 eV 
unless otherwise noted. It is convenient to fix the SDOS at some temperature T', and then solve Eq. 3 for 
with self-consistency requiring T{^th,T') = T'. 

When the OZ form of x is assumed, N_ becomes a constant, which we denote and Eq. 3 leads to long range 
order at T = 0 only, with correlation length ^th/d = 1 /given by 

ithqc = ( 11 ) 

where 5 = Uc — Usp^ Uc = l/xo(Qo: 0), T 2 = iiA 2 \/^Tuo?^ Ai is the coefficient of and Qc is a wave number cutoff. 

The origin of the anomalous region H for undoped LSCO, Eig. 2(b), can be readily understood from Eig. 2(a), 
where for T < 1500iF there is excess SDOS weight near Vc, leading to a strong peak (inset) as T ^ 0. This feature 
represents the development of the (tt, 7r)-plateau. Eor the small-t' materials, the susceptibility on the plateau remains 
parabolic, but with a small curvature and a sharp cutoff. This leads to an additional contribution to the SDOS of the 
form N_ = Np if X_ < X^. Then the integral of Eq. 3 becomes 

I = {Na- Nb5)ln{^ + 1 ) + NbXe + Npln{^ + 1 ). ( 12 ) 

where Xc is the maximum of X_. The last term, which we denote Ip, has two distinct limits. If Xp « S, then 
Ip NpXp/6. This is equivalent to assuming that the plateau has a flat top, in which case N- can be represented as 
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a ^-function, N_ = NpXpS{X_ — Uc), with NpXp the excess height of N at threshold, proportional to the area of the 
plateau at T = 0. As T decreases, S also decreases, and the opposite limit for Ip becomes appropriate, ^ << Xp, in 
which case Ip = Npln{Xp/5). From Eq. 3^ Ip ^ 1/5 translates into ^th r\j 1/T^/^, dotted line in Fig. 2(b). The actual 
anomaly roughly follows this line, but the correlation length growth is actually reversed. 

The present Gz — Wz — Tz model is the simplest which captures the essential physics of the pseudogap. For a 
more quantitative comparison with experiment, two additional problems must be solved. First, an extension to fully 
self-consistent G — W — T may be needed to capture the splitting of the saddle VHS peak into the VHS exciton and 
the residual continuum part. Secondly, the term in Eq. 1 proportional to must be included to describe quantum 
fluctuations. Since we are here primarily interested in the opening of the pseudogap at higher temperatures, we 
ignored it in our analysis. Finally, it will be interesting to explore if excitonic superfluidity exists, and how it is 
related to high-Tc superconductivity. 
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I. GENETIC ALGORITHM FOR CUPRATES 
A. A note on normalization conventions 

In this paper susceptibilities are calculated at in¬ 
termediate coupling level, using Z-corrected first- 
principles dispersions to describe the coherent part of 
the quasiparticle-GW dressed carriers.^ The results can 
be expressed in two different ways. If the susceptibility 
Xo is calculated based on the density functional theory 
(DFT) dispersion eo{k), then the phase diagrams are ex¬ 
pressed in terms of renormalized temperatures Tz = T jZ 
and Hubbard I/’s, Uz = ZU, where Z is the QPGW 
renormalization factor. This is used in the Supplemen¬ 
tary Material. For convenience in comparing with ex¬ 
periment, in the main text we calculated x based on 
the dressed dispersion ez = Zeo, assuming a doping- 
independent Z = 0.5. 

Also, when referring to Figures or Equations in the 
main text, we will use a leading ‘M’, as in Fig. Ml. 


B. Cuprate algorithm 

To understand differences between cuprates, and to 
gauge their proximity to the pure Hubbard model, it is 
necessary to add a third axis to the usual T- and doping- 
phase diagrams. Here we introduce a technique to reduce 
differences in band structure to the variation of a single 
hopping parameter. 

While DFT has been liberating in that it provides a 
parameter-free starting point for studies of materials, it 
can also be constraining, in that each material is defined 
as a single point in material parameter space, and it is 
hard to define connections between the members of a 
family of materials, such as cuprates or pnictides. What 
is needed is a form of materials genome, with an algo¬ 
rithm defining the relatedness of materials, the evolution 
from one material to another, and the probable route for 
extending the family to new materials. Here we propose 
such an algorithm for the cuprates - a ‘cuprate genome’. 
The cuprates are particularly simple, in that their prop¬ 
erties are dominated by a single Cu 02 band crossing the 
Fermi level, so their parameter space consists of hopping 
parameters, traditionally t, t', t", ..., which can be deter¬ 
mined, e.g., by a Wannier-type downfolding of the DFT 
bands. In fact, since t sets an energy scale, the band 
dispersion is determined by the ratios t'/t, ... In 



FIG. 1: Fluctuation fingerprints of cuprates. Maps of q- 
vectors of largest susceptibility for DFT-based model of LSCO 
(c) and Bi2201 (d) compared to t — t' — t" reference phase a 
[t" = 0] with t' = -0.23t (a) or -0.43t (b), and reference phase 
b [t" = with t' = -0.12t (e) or -0.258t (f). Recall, from 

Fig. Ml(c), that xq is the doping at T = 0, but we assume 

X ~ Xq. 


genetics, a family would be defined as a group lying on 
a low-dimensional surface in parameter space. The sim¬ 
plest genetics is associated with changes of a single gene, 
called Mendelian. Similarly, the simplest material family 
would lie along a single line in material space. We show 
that the cuprates correspond to such a Mendelian family, 
by introducing reference families depending only on t'/t, 
and show that the cuprates lie along a particular 
cut in this 2D-space. 

To do this, we introduce a definition of equivalence: 
ideally, two materials are equivalent if they have the same 
phase diagram. However, the real phase diagrams are 
defined by interactions beyond the current scope of DFT 
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FIG. 2: Mean-field phase diagrams of cuprates, for 

DFT-based models of LSCO (a) and Bi2201 (b). Black dots 
indicate mean field phase boundaries for different Uz values, 
including Uz = (down triangles), 1 (diamonds), or 0.5 eV 
(up triangles). 


theories. The best we can do is to say that they would 
have the same phase diagrams in the context of DFT. 
For weak and intermediately correlated materials, the 
ordered phases are generally determined from the Stoner 
criterion, and hence correlated with peaks in the Lind- 
hard susceptibility at cj = 0, which can readily be calcu¬ 
lated. 

Figure 1 is a map of the g-vectors associated with the 
largest susceptibilities, calculated for a variety of disper¬ 
sions. The left-hand frames (a), (c), (e), refer to LSCO, 
the right-hand frames (b), (d), (f), to Bi2201. The mid¬ 
dle row (c), (d), uses the full DFT dispersion, compared 
to two reference families, with (a) cuts defined by t" = 0 
(top row, (a), (b)) or (6) the Pavarini-Andersen cut^, 
t" = —t'/2 (bottom row, (e), (f)). We find that refer¬ 
ence family b semi-quantitatively reproduces the phase 
diagram for both cuprates. On the other hand, setting 
t" = 0 works well for LSCO, Fig. 1(a), but fails for Bi2201 
Fig. 1(b). We further find that a reasonable choice of 
t'/t is the one which matches the Van Hove singularity 
(VHS) doping of the full DFT dispersion. In these figures 
we fix an energy scale by taking t = 0.4195 eV, appropri¬ 
ate for LSCO^. We note that both reference families have 
the additional desirable feature that they evolve from the 
state with P = 0 - i.e., the original Hubbard model. 

This confirms our earlier finding that LSCO, Fig. 1(c), 
has a very different phase diagram from most other 
cuprates, here represented by Bi2201, Fig. l(d).^ In the 
phase labeled ‘ring’, the susceptibility is approximately 
constant along a ring of g-vectors surrounding (tt, tt). We 
note in particular that antinodal nesting (red-shaded re¬ 
gion) emerges naturally in most cuprates at higher dop¬ 
ing, but is absent in LSCO. Shown also in Fig. 1 are 
Tvhs (dark red dashed line) and the shift of the DOS- 
peak (thin green dashed line), discussed in the main text 
for LSCO. 

In Fig. 2, we replot the DFT fluctuation maps, adding 
the mean-field phase boundaries, based on a Stoner cri¬ 


terion, 

Uzxoiq,oj = 0) = ^- ( 1 ) 

In all cases, there is a critical value of U = Uq below 
which there is no phase transition at x = 0. For U < Uq, 
the phase boundary is dome-shaped, but with a peak 
shifted away from the T = 0 value of xvhs lo a value 
close to the line xvhs{T), red-dashed line. For U > Uq 
the phase boundary rapidly evolves to a strong coupling 
form where the highest T transition occurs near half fill¬ 
ing. This is consistent with our finding in the main text 
that there is an apparent shift of the VHS toward x = 0 
as T increases. In all cases, at sufficiently high T FS 
nesting is washed out and the dominant instability is at 
q = (7r,7r). 

An advantage of using the genetic algorithm is that it 
lets us systematically vary parameter space, to study the 
evolution from one cuprate to another, as in Fig. M2, 
or to extrapolate the parameter space to better under¬ 
stand the origin of various phases. Figure 3(a) shows 
the cuts in V — t"-space associated with reference fami¬ 
lies a {t" = 0) [red line] and b {t" = —t'l2) [blue line], 
where the light-green shaded segments show the range 
over which ^ is significantly depressed. Fig. M2. The 
dotted lines indicate ground state transitions at x = 0, 
from (tt, tt) to (tt, tt — S) to {6, 6) order. The densities- 
of-states (DOS) shown for the two families in Figs. 3(b) 
and 3(c) provide insight into the origins of the crossovers. 
For reference family a the DOS reveals two crossovers as¬ 
sociated with a change from 2D to ID nesting - note that 
the DOS at t' = —1/2 shows the characteristic 1/x/IE) 
divergence expected for a ID system. Near t' = —0.5t, 
the DOS has two peaks, one at the saddle-point VHS 
found in 2D materials, the other at the leading edge, as 
found in ID systems. For t' < —0.46t, the leading edge 
is larger, suggestive of weakly-coupled chains, while for 
t' > —0.2St the leading edge peak turns into a local min¬ 
imum. This quasi-ID behavior is suggestive of the ne¬ 
matic phase found in cuprates. Remarkably, t' = — 0.28t 
approximately coincides with the C-I transition found in 
Fig. M2 to lead to a severe ^ depression. This explains 
the extended T-range of ANN [((^, (^)] nesting found in 
Fig. 1(b), which is not characteristic of cuprates. 

The origin of the ID-behavior can be easily under¬ 
stood. For reference family a, the electronic dispersion 
has the form E = —2t{cx + Cy) — At'cxCy —2t when 
t' = —0.5t and either kx or ky = 0. [Here q = cos{kia), 
i = x^y.] That is, at the VHS the FS reduces to just 
the X and y axes. This explains the strong role of t", 
which breaks up this perfect nesting, and shifts the dop¬ 
ing of optimal ANN nesting away from the VHS-doping. 
This causes the FS to split into two pockets, with the 
VHS corresponding to the topological transition where 
the pockets split, inset to Fig. 3(b). The onset of this 
transition lies close to the doping for Bi2201 (light blue 
DOS cur). 
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FIG. 3: Reference families of the cuprates, (a) Refer¬ 
ence cuts in t' — t"-space. Dotted lines indicate ground state 
transitions at x = 0. Light green line segments indicate the 
ranges over which ^ is significantly depressed, Fig. M2, with 
maximal depression indicated by yellow diamonds. Violet tri¬ 
angles [green squares] = reference states for LSCO [Bi2201]; 
colored circles indicate t'-yalues for the corresponding DOSs 
in frame (b) [upper line] or (c) [lower line], (b) DOS for 
several values of t' for reference family b (t" = —1'/2). As 
VHS moves from right to left, t'/t — 0 (red curve), -0.1 (or¬ 
ange), -0.2 (yellow-green), -0.25 (green), -0.258 (light blue), 
-0.3 (blue), -0.4 (violet), and -0.5 (black). Inset: correspond¬ 
ing Fermi surfaces at VHS. (c) DOS for several values of t' for 
reference family a {t" = 0). As VHS moves from right to left, 
t'/t — 0 (red curve), -0.3 (orange), -0.4 (yellow-green), -0.45 
(green), -0.47 (light blue), -0.48 (blue), -0.49 (violet), and -0.5 
(black). Inset: leading edge of DOS, near band bottom for 
(from bottom to top) t'/t — -0.22, -0.27, -0.28, -0.29, -0.30, in 
regime where it crosses over from a local minimum to a local 
maximum. 


II. COMPETITION BETWEEN FS AND VHS 
NESTING 


ity that enters into the Stoner criterion, 
^ Ck- eu+q 


( 2 ) 


The denominator is largest if e/e = e/e+g, but in this case 
the numerator vanishes unless is at the Fermi level. 
Even in this case, the term contributes little weight to 
the susceptibility integral unless the g-shifted FSs are 
tangent, which is generally satisfied if g = 2 /ci?, giving 
rise to a folded map of the FS - the ridge of FS nesting. 

This can be seen by comparing a map of Xq, Fig. 4(a), 
with the corresponding map of Xq, Fig. 4(b). The plot is 
appropriate for a related cuprate, Bi 2 Sr 2 Cu 06 (Bi2201), 
at X = 0.13 and cj = 5 meV, showing in the process 
that VH nesting also produces a significant effect in 
the Bi-cuprates. Due to the energy (i-function, Xo con¬ 
tains mainly the FS contribution^. [The extra feature 
at F is the DOS contribution.] One sees a very simi¬ 
lar ‘inverted- 7 ’ shaped ridge in both Xo and Xo? which 
is a g = 2/ci? map of the Fermi surface^. In Fig. 4(d), 
we plot Xo (i*cd line) and Xo (blue line) along a cut 
F ^ (tt, 0) ^ (tt, tt) ^ F in momentum space. It can 
be seen that, away from F, the peaks in both compo¬ 
nents fall at the same g-values, and the relative weights 
of the peaks is comparable. Thus, the anti-nodal peaks 
near F are stronger than the near-( 7 r, tt) peaks. However, 
Xo is largest on the near-( 7 r, 7 r) peaks, since the FS con¬ 
tribution is riding on a large peak associated with the 
bulk contribution to the susceptibility. 

This bulk contribution arises as follows. Even when 
e/e 7 ^ c/e+g, there will be a finite contribution to the sus¬ 
ceptibility, reduced by the factor l/(e/e — e/^+g), as long 
as the two energies lie on opposite sides of the Eermi 
level (for T = 0). This contribution will be especially 
significant near a VHS, where there is a large local DOS 
(LDOS). Indeed, this contribution will remain large even 
as X is doped away from xvhSi since while the LDOS on 
one side oiEp decreases, that on the other side increases, 
so that the product remains little changed until the VHS 
peak is quite far from Ep- This bulk part provides a 
smoothly varying background, with hills that can shift 
the balance of the ES nesting, and in special cases can 
lead to commensurate nesting away from the ES nest¬ 
ing vector. Here we try to quantify the dominant role of 
VHS nesting in creating this background, thereby pro¬ 
viding a deeper understanding of the VHS vs ES nesting 
competition. 


A. Origin of FS and VHS nesting 

Here we attempt to separate the two contributions to 
the susceptibility, from bulk and Fermi surface, and to 
clarify the origin of the bulk contribution. To our knowl¬ 
edge there is no standard way to do this, so we have 
developed simple schemes. 

Consider the expression for the 00 = 0 bare susceptibil¬ 


B. Separating FS and VHS nesting contributions 

Two approaches are taken to eliminate the FS nest¬ 
ing features to see the background susceptibility by it¬ 
self. First, raising T (violet long-dashed line in Fig. 4(d)) 
washes out the sharp FS features first, while leaving the 
background features behind. Secondly, in frame (c), we 
plot the difference between Xo and suitably scaled quan¬ 
tity proportional to Xm this is also plotted as the green 
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FIG. 4: Extracting the bulk xo contribution, (a) Sus¬ 
ceptibility Xo of Bi2201, X = 0.13, plotted in first Brillouin 
zone, at T = 0, cc = 5 meV. Color scale from minimum (blue) 
to maximum (white), (b) As in (a), but for xi- (o) Difference 
plot of (a) and scaled version of (b). (d) Cuts of suscepti¬ 
bility from earlier frames plotted along high symmetry lines, 
with red curve representing data in frame (a), blue curve for 
frame (b), and green curve for frame (c). Shown also is the 
curve corresponding to frame a, but at T = lOOOK (violet 
long-dashed line) scaled by 1.4 and shifted by -0.045 eV“^ to 
match other curves near (7r,7r).. 


FIG. 5: Coherent-incoherent crossover at x=0: role 
of VHS nesting. Deconvolution of the susceptibility for 
/c-states in the (tt, 0)-patch at several temperatures T. Con¬ 
tribution to the susceptibility of Xo, 2 (^,e/c), Eq. 3, is plot¬ 
ted along ^-cuts near (tt, tt) with the color scale ranging from 
white for a large contribution to blue for negligible contribu¬ 
tion for: 0 (a), 25K (b), 50K (c), 250K (d), 500K (e), and 

2000K (f). 


curve in frame (d). Since x'o Xo Kramers-Kronig 
transforms, they have different lineshapes, but since the 
FS map is nearly one-dimensional, a large part of it is still 
cancelled in taking the difference. From frame (d) it can 
be seen that both approaches lead to similar results for 
the bulk susceptibility, and that this contribution pro¬ 
vides the dominant contribution to virtually the whole 
(tt, tt) plateau structure. 


C. Susceptibility Analysis 

Finally, we demonstrate the connection between the 
bulk contribution to xo and the (tt, tt) bosonic-VHS (b- 
VHS), by deconvolving the bare susceptibility, Eq. 2, into 
its various components. To understand this procedure, 
recall that the VHS in the electronic dispersion arises 
when the Fermi surface crosses the Brillouin zone bound¬ 
ary near (tt, 0) and related points, crossing from hole-like 
to electron-like. The electron-hole pairs associated with 
the VHS come in two kinds: intra-VHS pairs, where both 
carriers lie near the same VHS, so their relative momen¬ 
tum is near q = (0,0) = F, and inter-VHS pairs, where 


an electron near (tt, 0) and a hole near (0, tt) have rel¬ 
ative momentum near (tt, —tt). The latter pairs carry 
more weight, forming a b-VHS near q = (7r,7r). Thus, to 
confirm that the susceptibility peak at (tt, tt) is assiciated 
with the b-VHS, we analyze Eq. 2 for xo, assuming that 
the electron is associated with and the hole near 6/.+^. 
We must demonstrate that, when q ^ (7r,7r), the domi¬ 
nant contribution to xo has k near the VHS at (tt, 0) [or 
(0,7r)] and ek near Evhs- 

We divide the Brillouin zone into four patches, around 
the F, (7r,0), (0,7r), and (7r,7r) points, and calculate sep¬ 
arately the susceptibility due to each patch. In LSCO, 
the plateau is predominantly composed of equal contri¬ 
butions from the (0, tt) and (tt, 0) patches, as expected 
for a VHS. We then explore the distribution of values 
of e/e which contribute to the susceptibility in the (tt, 0) 
patch at X = 0, Fig. 5, rewriting Eq. 2 as 

xo{q) = EE Xo,i(g,efe), (3) 

i=l,4 

where i runs over the four quadrants F, (tt, 0), (0, tt), and 
(7r,7r). 

There is a clear crossover from a dominant near-FS 
contribution (e/e = 0) at low T to a dominant contri- 
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bution at high-T coming from the VHS, caused by the 
breakdown of Pauli blocking, Fig. Ml(b). At high T, 
the dominant contribution lies along a flat line centered 
at energy proportional to ksTyns = Ep — Eyns- The 
same crossover is evident in the susceptibility, as the co¬ 
herent, FS-related features present at low T are washed 
out, leaving only a featureless (incoherent) (tt, tt) plateau 
at high T, Fig. 6 . As doping increases, the same pattern 
is repeated, but with the flat line moving toward Ep as 
Tvhs decreases, and the crossover from coherent to inco¬ 
herent shifts to lower T, consistent with the pink shaded 
region in Fig. Ml(c). The boundaries of this shaded re¬ 
gion are Ty 7 / 5/2 and Ty 7 / 5 / 3 , where Typs is plotted as 
the dark red dashed curve in Fig. Ml(c). 

While this explains what happens at T > Typs^ 
question remains why the broad peak at (tt, tt) persists 
to lower T [e.g.. Figs. l(a,c,e) at x=0]. To understand 
this we first look at what happens when Pauli block¬ 
ing is ostensibly absent - when Ep = Eyps, Fig. 6 (b). 
While the VHS peak is clearly present at T = 0 (vi¬ 
olet line), finite temperature does not lead to simple 
thermal broadening but actually to a split peak. This 
is a thermally-induced Pauli blocking effect: for |ek|, 
kk+ql < T, /k = /k+q = 1/2 and A/ = 0. Since the 
(tt, tt) susceptibility is pinned at cj = 0, the Pauli block¬ 
ing turns on first at (tt, tt) and then spreads at higher T. 
This has a number of consequences: (1) As x shifts to¬ 
ward half filling, the peak splitting turns on at higher Ts, 
Fig. 6 (a) - this is the lower branch of the (tt.tt) — (tt, tt — (5) 
crossover, plotted in Fig. Ml(c). (2) This is the origin of 
the ring phase, discussed above in relation to Figs. l(d,f), 
where the peak susceptibility shifts from (tt, tt) to a ring 
of states encircling ( 7 r, 7 r). (3) Finally, since the integral 
of / is fixed by the net doping, when a thermal blocking 
turns on, some other states will be unblocked at nearby 
energies. Thus, the persistence of the net (tt, tt) peak at 
T < Typs still reflects the underlying presence of the 
VHS, but is more and more controlled by states away 
from ( 7 r, 7 r) and closer to Ep. This evolution can be 
clearly traced in Fig. 5. It is in this sense that we can 
refer to this VHS as hidden. 

The low-T fadeout of the b-VHS is controlled by the 
Pauli blocking factor. 


^fk,Q 


sinhx- 

coshx- -h coshx+ ’ 


( 4 ) 


where X- = eQ-{k)/kBT, Xp = (eQ,+ (/c) - Ep)lkBT. 
When the electronic VHS is at the Fermi level, Xp = 
0 and Afk^Q = tanhx_/ 2 , while at lower doping 
the b-VHS is exponentially suppressed, coshx+ ^ 
exp{TypsIT)I2^ which equals 1 when T ^ Typs/ln{2). 



FIG. 6: Lindhard susceptibility of LSCO near the 
VHS doping. Xo(^,Ct; = 0) is plotted for a series of tem¬ 
peratures (from top to bottom) T = 0 (violet), 10 (blue - 
frame (b) ouly), 25, 50, 100, 200, 400, 1000 (greeu), 2000, 
aud 4000K at dopiug x = 0.18 (a) aud 0.20 = xyps (b). 


III. INTERPRETATION OF 
SELF-CONSISTENCY 

A. Relation to Bose-Einstein Condensation 

In the conventional picture of Bose-Einstein conden¬ 
sation (BEC), the total number of bosonic modes must 
satisfy 

E = J dujnBE{uj)DB{uj), (5) 

where ube the BE distribution function and Db is the 
bosonic density of states. As T decreases, the bosonic 
chemical potential p must adjust to maintain the equal¬ 
ity. However, the largest fi can be is zero, for which 
^be the Planck distribution, and the RHS of 

Eq. 5 takes a fixed value. Eor lower T the equation 
can only be satisfied by putting a macroscopic number of 
bosons into the lowest energy mode. 

In the present situation, Db = ^qX"^ ^p can be 
approximated as ^ Tjuj, in which case Eq. 5 becomes 

N ( dw ^ = 0 )- ( 6 ) 

q q 

Comparing this to Eq. 5 in the main text, self-consistency 
has produced a number N^ff = NttX/Aq of excitons 
which can subsequently Bose condense. Whereas for 3D 
BEC the macroscopic occupation of the lowest bosonic 
mode must be added by hand, in 2D the self-consistent 
condition automatically causes the occupation of that 
mode to diverge as T ^0. 
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B. Entropy and Dissipation 

As originally noted by McMillan, anomalously large ra¬ 
tios of 2A/kBTc for any phase transition are generically 
associated with a dominance of bosonic entropy at finite 
T, associated with low-lying bosonic fluctuations. But 
from the fluctuation-dissipation theorem, the net fluctu¬ 
ations at a given frequency uj are equal to the dissipation 
D{uj) at that frequency, 

-D(w) = '^x" {q,u:){nBE{(^) + 1/2) 
q 

(7) 

q 

so the right-hand side of Eq. 6 is just the total dissipa¬ 
tion at temperature T. In an ordered phase near T = 0, 
the dominant fluctuations will be associated with Gold- 
stone modes, but to show this in a calculation requires 
extending the present results to full self-consistency be¬ 
tween self energy and vertex corrections, as in a parquet 
calculation. 


IV. MODE COUPLING IN BI2201 

In Bi2201 and most other cuprates, a larger ratio \t'/t\ 
expands the (tt, tt)- plateau, leading to a rich phase dia¬ 
gram, Fig. 2(b). In this case the b-VHS is dominant only 
at high T, while FS-related susceptibility features cause 
a crossover of the (tt, tt) susceptibility from a maximum 
at high T to a local minimum at lower T. This leads to 
a different kind of strong mode coupling at intermediate 
T, best exemplified by the ring phase. Figs. l(c,d). In 
this phase, the leading-edge SDOS has a ‘Mexican hat’- 
like divergence [Main Section II.C], Fig. 7(a) near 2000K, 
leading to a slow growth in ^ ^ Fig. 7(c). Near the 
onset of ring order, Tz = 32507^, the correlation length 
has a nonmonotonic T-dependence, similar to that seen 
in LSCO, Fig. M2(b). 

While many pseudogap-like features are similar to 
LSCO, Bi2201 has some novel properties associated with 
the ANN instability, which becomes a CDW when charge 
order is included in the model. Since the ANN suscepti¬ 
bility peak has a much larger intrinsic curvature, mode 
coupling effects are weak. Hence the ANN susceptibility 
peak grows rapidly at low temperature, and as doping 
increases, it can surpass the (tt, 7r)-peak, thereby sharply 
cutting off the growth of near-(7r,7r) fluctuations. This 
is illustrated in Fig. 7 for x = 0.20. Notably, since this 
peak is not associated with strong mode coupling it has a 
weak SDOS with a very strong T-dependence, Figs. 7(b) 
and 7(d), crossing the NNN peak near Tz = 200K. The 
large subsequent growth of the ANN correlation length 
would ensure the rapid establishment of a true long-range 
3D ANN order. However, corresponding to the large T- 
dependence of Xann there is a similar strong disorder 
dependence - particularly in the presence of short-range 



FIG. 7: Structures in Bi2201 correlation lengths. Tem¬ 
perature dependence of: (a) SDOS N- for t — t' — t" reference 
state of Bi2201, x = 0.20. Data sets are offset for clarity; 
from highest to lowest, Tz = 4000 (red), 2000 (violet), 1000 
(blue), 750 (green), 500 (brown), and 200K (light blue), (b) 
Continuation of SDOS in (a) to lower Tz, showing evolution 
of leading edge. From top to bottom: Tz = 200 (light blue), 
100 (red), 50 (blue), and lOK (green). For each curve, a 
partial SDOS near the ANN peak is superposed to show the 
SDOS associated with the ANN peak whose onset is marked 
by arrows, (c) Correlation length vs Tz- Blue dotted line 
represents ^ ~ T~^ . (d) Susceptibility of x{q) peak vs Tz for 
q near the (tt, 7r)-plateau (red line) or for ANN nesting (blue 
line). 


near-(7r,7r) order, which should cause the ANN corre¬ 
lation length to saturate. Nevertheless, the transition 
should be manifest as a coherent-incoherent crossover in 
the overdoped regime, which shifts to higher T as dop¬ 
ing increases. Remarkably, just such a feature has been 
found in ARPES studies.^ Combining this crossover with 
the T* crossover at lower doping produces the V-shaped 
boundary of the ‘strange metal’ phase associated with 
linear-in-T resistivity. 


V. VHS EXCITONIC INSTABILITY AND 
KONDO-LIKE PHYSICS 


The physics of semimetals is strongly influenced by the 
simultaneous presence of both holes and electrons, and 
the attractive interaction between them. In lightly doped 
semimetals this can lead to exciton condensation, an ex¬ 
tension of the excitonic instability found in small gap 
semiconductors^’^. [See Ref. 9 for a recent experimental 
example.] While screening modifies excitons in semimet¬ 
als and doped semiconductors, excitonic effects do not 
disappear, but shift from ^-function bound excitons to 
excitonic resonances with a finite width, which persist to 
high dopings and show strong, even resonant interaction 
with phonons^^. Indeed, Mahan has proposed that ex¬ 
citon formation in a doped semiconductor is resonantly 
enhanced by the presence of a sharp Fermi surface, in 
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analogy with the edge singularity in x-ray spectrad^ 

The strong effect of excitons in doped semiconductors 
raises new questions about the role of VHSs in metals. 
At a saddle-point VHS in a two-dimensional (2D) metal, 
or an Mi VHS in a 3D metaff^, the topology change of 
the Fermi surface (FS) leads to a large coexisting popu¬ 
lation of electrons and holes. Indeed, resonant interhand 
excitons associated with the Mi VHS (VH excitons) are 
known to dominate the optical spectra of metals. Sim¬ 
ilarly, in cuprates, an intrahand Van Hove exciton can 
arise near an electronic saddle-point VHS^^’^^. Just as in 
conventional exciton Bose condensation, the mean-field 
transition corresponds to the T where [quasi] excitons 
form, but since these are nearly localized, there is a large 
spread in g, leading to incoherent VHS nesting. The real 
transition lies at much lower T, when all pairs condense 
into the lowest energy g-state. This excitonic VHS con¬ 
tinues to play a role even when the electronic VHS is 
away from leading to a peak in x'o when either T or 
(j is Ep — EvhS‘ This is due to scattering between 
states near the VHS and states near Ep^ and is similar 
to Kondo physics, with the VHS states playing the role 
of /-electrons. 

To clarify the connection to Kondo physics, we write 
a variational ground state for the singlet exciton as 

I'lpG >= + Vfe4+Q,<TCfc,<T) I-FCtt.tt) >, (8) 

k,a 

where ^ [ck,a] are creation [annihilation] operators for 
an electronic state of momentum k and spin <j, and 
Uk^ Vk are variational parameters. Here > is a 

filled Fermi sea of the first (tt, 7r)-superlattice Brillouin 
zone, and the excitonic order parameter is Acdw =< 
^k+Q T ^I+Q T c.c. >. This model is a direct 
generalization of Yosida’s variational calculation of the 
screening of a Kondo impurity to a single-band Kondo 
lattice model.Moreover, to form a localized exci¬ 
ton wave-packet, the restriction g = Q in Eq. 8 must 
be dropped, and many g-vectors summed over, account¬ 
ing for the bosonic entropy of the model - indeed the 
SDOS quantifies the number of near-Q states that con¬ 
tribute to the entropy. Since the Fermi velocity vp = 
at a VHS, these excitons are nearly localized. This rein¬ 
forces the connection of the coherent-incoherent crossover 
to Kondo physics^^: in Main Fig. M2(b), region I corre¬ 
sponds to the buildup of a Kondo lattice of nearly lo¬ 
calized spins, while region H is a regime of crossover to 
itinerant physics. 

VI. RELATION TO OTHER CALCULATIONS 

A comparison between our QPGW model and other 
‘DFT-h’ calculations was given in Ref. 1; here we dis¬ 


cuss this comparison specifically in regard to pseudogap 
physics. A number of groups are working on MBPT- 
type extensions of DFT calculations. However, DMFT 
and its cluster extensions are limited by a very low res¬ 
olution in momentum g. In particular, this means that 
even if they capture Mermin-Wagner physics, they are 
typically limited to ^ [Fig. 50 of Ref. 1]. Fur¬ 

thermore, in averaging over large patches in momentum 
space, they lose the ability to resolve ordinary nesting 
instabilities, such as the ANN instability responsible for 
cuprate CDWs [Fig. 49 of Ref.l]. On the other hand, 
the fact that they can see the pseudogap, exen when av¬ 
eraging over a quarter of the Brillouin zone is consistent 
with our finding that the pseudogap is spread over many 
separate modes in g. 

However, some DMFT extensions [dynamic vertex 
approximation^^^^^, dual fermion^^’^^, and one-particle 
irreducible approaches'^] seek to calculate the two- 
particle vertex, which allows them to calculate the 
momentum-dependent self energy, and these approaches 
are able to reproduce Mermin-Wagner physics and vertex 
corrections. We note in particular that a recent dynam¬ 
ical cluster approximation (DCA)^^ calculation sorted 
out the contribution of various channels of two-particle 
scattering to the self energy (‘fiuctuation diagnostics’), 
and found that spin fiuctuations are the origin of the 
pseudogap, while charge and pairing fiuctuations play a 
marginal role, consistent with our results. This has also 
been found experimentally^^. 

We further note that the parquet equations^^’^^ could 
be quite useful, if the recently discovered divergences^^ 
can be overcome. The mutual self-consistency between 
the self-energy and vertex corrections will be necessary 
to extend the current results to low temperatures. 

Finally, in many strong-coupling calculations, it 
is implicitly assumed that Fermi surface nesting is 
irrelevant^^, and there is a special role for commensu¬ 
rate (7r,7r) magnetic order. Thus, in the Yang-Rice- 
Zhang phenomenological modefi^, the nesting is tied to 
the magnetic Brillouin zone boundary. The commensu¬ 
rate phase we have found provides a clear DFT-based 
explanation for that effect, as well as a clear picture of 
how this evolved into a more Fermi-liquid-like phase at 
higher doping or larger t'. 

After our paper was posted [arXiv: 1505.04770], we be¬ 
came aware of another paper^^ that also finds that strong 
bosonic fiuctuations over a ‘quasidegenerate distribution 
of 2pF ordering wave vectors’ drives the pseudogap, al¬ 
though they find it is predominantly in the charge sector, 
forming a ‘resonant Peierls excitonic state’. 
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